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EFFECT  OF  TRANSDUCER  SIZE  ON  THE  STATISTICAL  PROPERTIES  OF 
A  SIMULATED  TURBULENT  WALL  PRESSURE  FIELD 


By 

Robert  L.  Ash^  and  Mehdi  R.  Khorrami2 
\\  SUMMARY 

The  root  mean  square  pressure  was  calculated  with  synthetically 
generated  pressure  signals  for  transducers  of  different  size.  The 
transducer  diameters  were  in  the  range  of  0.02  <  d  <1.18  at  a  Reynolds 

number^Ra^y  of  12,250.  The  lower  limit  of  the  diameter  simulated  here 
was  much  smaller  than  diameters  reported  previously  by  other  workers. 
Calculated  root  mean  square  pressure  levels  developed  in  the  present 
work  were  in  good  agreement  with  the  measured  experimental  data. 

In  order  to  check  the  simulated  pressure  signals,  the  power  spec¬ 
trum  and  two  point  correlation  were  calculated  as  well.  The  power  spec¬ 
tra  obtained  with  the  smaller  transducers  were  higher  in  energy  level  at 
high  frequencies.  The  shape  of  the  power  spectra  were  similar  and  over¬ 
all  agreement  with  experimental  data  was  reasonable.  The  two  point  cor¬ 
relations  were  obtained  for  separation  distances  of  .66,  and 

13.37.  The  same  convection  and  decay  effects  as  the  experimental  data 

/ 

were  observed  with  the  calculated  correlations.  Again,  overall 
agreement  with  experimental  result  was  quite  good. 
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Old  Dominion  University,  Norfolk,  Virginia  23508. 

2  Graduate  Research  Assistant,  Department  of  Mechanical  Engineering  and 
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I.  INTRODUCTION 


Turbulent:  wall  pressure  data  are  required  for  Che  analysis  of  engi¬ 
neering  problems  related  to  noise  production  on  aerodynamic  bodies,  as 
well  as  vibration  and  panel  flutter  studies  (Dowell  [1])*.  In  recent 
years,  the  identification  of  compliant  walls  which  are  receptive  to 
turbulent  pressure  fluctuations  has  been  of  particular  interest  (Bush- 
nell  et  al.  [2]  and  Buckingham  et  al.  [3]).  It  has  been  speculated  that 
a  compliant  surface  might  be  capable  of  interacting  with  a  turbulent 
flow  to  produce  a  net  drag  reduction.  When  local  instantaneous  pressure 
data  are  required,  the  resolution  capability  of  even  the  smallest 
pressure  transducers  becomes  a  serious  limitation  because  they  average 
the  pressure  field  spatially  (Corcos  [4]  and  Emmerling,  Meir  and 
Dinkelacker  [5]).  The  purpose  of  this  investigation  has  been  to  utilize 
a  synthetic  pressure  field  to  study  the  effect  of  transducer  size. 

Instantaneous  wall  pressure  fluctuations  can  be  generated  either  by 
employing  large  arrays  of  pressure  transducers  with  prescribed  spatial 
resolution  in  an  experiment  or,  by  solving  the  full  three-dimensional 
Navier-Stokes  equations.  Although,  both  methods  can  be  more  accurate 
than  a  simulation,  they  represent  a  formidable  task  in  terras  of  time  and 
expense  and  are  restricted  to  a  particular  temporal  and  spatial  resolu¬ 
tion.  The  Honte  Carlo  technique  employed  in  the  present  work  is  capable 
of  resolving  arbitrarily  in  time  and  space. 

The  contributions  from  small  scale  motions  (with  length  scales  on 


*  The  numbers  in  brackets  indicate  references. 


the  order  of  —  ,  where  UT  is  the  friction  velocity)  to  the  wall  pres¬ 
sure  cannot  be  measured  currently  because  no  pressure  transducers  are 
sufficiently  small.  However,  during  the  late  sixties  and  early  seven¬ 
ties  the  importance  of  high  frequency  pressure  fluctuations  on  the 
measurements  of  root-mean  square  pressure  (P  )  was  established  by 

Blake  [6]  and  Emmerling,  Meir  and  Dinkelacker  [5]  and  a  rough  estimate 
of  their  contribution  to  P„  u  „  was  given  by  Emmerling,  Meier  and 
Dinkelacker  [5J.  The  understanding  of  the  above  features  were  made 
possible  by  employing  pressure  transducers  vrtiich  were  miniature  when 
compared  to  earlier  experiments.  Data  generated  from  these  sensors 
indicated  that  the  R.M.S  pressure  increased  asymptotically  with  decreas¬ 
ing  transducer  size.  The  present  work  has  attempted  to  clarify  the 
above  suggestion  by  using  synthetic  pressure  fluctuations  to  study  the 
effect  of  a  controlled,  repeatable  pressure  history  on  the  measure¬ 
ments  from  different  sensor  sizes.  In  order  to  validate  the  synthetic¬ 
ally  generated  pressure  fluctuations,  other  important  features  of  a 
turbulent  boundary  layer  (such  as  power- spectrum  and  two  point  corre¬ 
lation)  have  been  calculated  as  well. 

In  recent  years,  a  unified  view  has  evolved  concerning  the  exis- 
tance  of  two  different  scales  of  disturbances  in  a  turbulent  boundary 

layer.  Some  disturbances  have  scales  which  are  comparable  to  displace- 
★ 

ment  thickness  (6  ,  which  can  be  related  directly  to  boundary  layer 

thickness)  and  other  disturbances  have  much  smaller  characteristic 

v 

lengths  which  are  on  the  order  of  — —  .  It  is  now  obvious  from  the  log- 

Ux 

log  plot  of  power  spectrum,  that  any  shift  in  the  high  frequency  region 
due  to  the  measurements  of  small  scale  eddies,  can  result  in  an  increase 


2 


in  R.M.S.  pressure  (Blake  [6]),  since,  PD  u  e  is  related  to  the  power- 

K« rl*  a» 

spectral  density  by: 

P2  *  2 /°°$  (a)  )doi . 

o 

Kraichnan  (8)  first  tried  to  calculate  important  features  of  a 
pressure  field  by  integrating  Poisson's  equation 

32U.U. 

v2p  =  "  p  Sx.Sx"!  » 

1  J 

over  a  finite  region.  In  a  subsequent  paper,  Kraichnan  [9]  assuned 
local  isotropy  in  planes  parallel  to  the  surface  and  used  a  crude  model 
for  mean  velocity,  to  estimate  the  contribution  from  turbulent  mean 
shear  interactions  (which  he  concluded  were  the  dominant  terms) .  He 
also  assumed  that  any  column  of  fluid  with  a  length  scale  comparable  to 
the  boundary  layer  thickness  was  statistically  independent  of  the 
adjacent  columns.  The  estimated  R.M.S  pressure  was  found  to  be  about 
six  times  the  wall  shear  stress.  Recent  experiments  by  Emmerling,  Meier 
and  Dinkelacker  (5]  and  Bull  and  Thomas  [10]  have  shown  that 
Kraichnan' s  PD  M  e  estimate  was  in  good  agreement  with  their  measure- 
considering  the  degree  of  approximation  he  had  employed. 

Corcos  [4]  integrated  a  corrected  power- spec trim  (his  method  will 
be  discussed  later)  and  was  able  only  to  give  the  ratio 


where  P2  is  the  integrated  power  spectrun  of  Willraarth  and  Wooldridge 


3 


[14].  The  integration  was  up  to  radian  frequencies  given  by  — —  ■  14.3, 

00 

which  was  the  highest  measured  frequency  by  Willmarth  and  Wboldridge 
[14].  He  indicated  that  the  contribution  from  higher  frequency  fluctua¬ 
tions  might  be  significant.  Emmerling,  Meier  and  Dinkelacker  [5]  have 
employed  an  optical  method  to  evaluate  the  fluctuation  field  at  the 
wall.  For  the  sake  of  comparison,  they  measured  the  R.M.S.  pressure 
with  microphones  mounted  flush  with  the  wall.  They  plotted  the  varia¬ 
tion  of  non-dimensionalized  P  with  non-dimensionalized  transducer 

R.M.S • 

diameter  and  compared  them  with  measurements  of  several  other  experi¬ 
ments.  Their  variables  were  made  dimensionless  with  respect  to  free 

v 

stream  dynamic  pressure,  q  ,  and  inner  variables  — —  ,  respectively. 

uT 

Their  plot  showed  a  dramatic  increase  in  R.M.S.  level  as  the  non-dimen¬ 
sionalized  diameter  was  decreased  below  p  150.  Emmerling,  Meier  and 

v 

Dinkelacker  [5]  have  indicated  that  when  the  scaling  of  transducer 

diameter  was  with  respect  to  the  displacement  thickness,  6*,  the  data 

did  not  collapse  onto  a  single  curve.  This  suggested  that  small  scale 

eddies  were  major  contributors  to  the  P  „  c  (or  power  spectrum)  of  a 

K.H.  b. 

turbulent  boundary  layer. 

In  order  to  understand  the  results  given  by  Emmerling,  Meier  and 
Dinkelacker  [5],  Bull  and  Thomas  [10]  conducted  an  experiment  where 
direct  power-sRectral  measurements  were  made  using  pinhole  microphones 
and  piezoelectric  transducers  with  the  same  diameter.  With  a  similar 
plot  to  that  of  Emmerling,  Meir  and  Dinkelacker  [5],  they  have  shown 
that  although  there  is  an  increase  in  P.  M  .  with  decreasing  diameter, 


4 


the  increase  is  not  as  significant  as  reported  earlier.  Bull  and  Thomas 

[10]  have  attributed  the  high  values  of  P  measured  by  Blake  (6) 

R.M.S. 

and  Emmerling,  Meier  and  Dinkelacker  [5]  to  the  local  flow  disturbances 

created  by  the  pinhole-microphones. 

Recently,  Bull  and  Langeheineken  [11]  reproduced  the  experiment  of 

Bull  and  Thomas  [10].  They  reinterpreted  the  transducer  size  parameter 
+ 

d  ,  which  was  Reynolds  number  dependent  and  given  by: 


d 


+ 


U 

ci  _x_  Re. 

L  Uco 


where  L  is  the  characteristic  length  of  the  flow  field  (pipe  radius  or 
boundary  layer  thickness) .  Their  plots  of^P^/q  vs  d+  and  Re  show 
similar  trends  for  condensor  microphones  and  piezoelectric  transducers. 
Therefore,  they  concluded  that  condenser  microphones  and  piezoelectric 
transducers  give  essentially  identical  results  (which  is  contrary  to  the 
explanation  of  Bull  and  Thomas  [10]).  However,  the  above  conclusions 
have  produced  other  questions.  Firstly,  since  the  transducer  diameters 
were  not  equal,  direct  comparison  with  Bull  and  Thomas'  [10]  results 
were  not  possible.  Secondly,  their  experiment  was  conducted  in  a  pipe 
flow  which  raises  questions  concerning  how  different  flow  fields  affect 
the  turbulent  wall  pressure  field.  Specifically,  the  question  of  whether 
the  pressure  gradient  imposed  in  a  pipe  flow  alters  the  statistical 
properties  of  the  pressure  fluctuations  must  be  addressed. 
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Most  of  the  semi-empirical  work  concerned  with  the  correction  for 
transducer  size  has  followed  the  method  first  outlined  by  Uberoi  and 
Kovasznay  [12].  They  developed  a  mapping  procedure  based  on  a  knowledge 
of  the  instrument  response  which  enabled  them  to  form  a  response  kernel 
for  the  measuring  system.  The  original  signal  (or  the  instantaneous 
signal,  since  they  have  assumed  the  instrument  response  is  perfect) 
could  be  retained.  But,  the  limitation  of  their  method  was  due  primari¬ 
ly  to  the  assumption  of  a  total  random  field  which  was  isotropic  and 
homogeneous.  The  established  fact  that  turbulence  has  a  "fading  memory" 
implies  that  it  is  neither  homogeneous  nor  isotropic.  Corcos  [4]  util¬ 
ized  the  work  of  Uberoi  and  Kovasznay  [12]  to  develop  a  similarity  model 
for  cross-spectral  density  in  the  form: 


r  0,5  ,n) 


where  a  / ai5  \  and  3/  torj_ \  are  similarity  functions  (based  on  the  ex¬ 


perimental  data  of  Willmarth  and  Wooldridge  [14])  for  the  longitudinal 
and  lateral  directions,  respectively.  He  assumed  a  uniform  response 
over  the  transducer  surface  which  had  a  kernel  function  equal  to  the 
inverse  of  its  area.  That  form  is  consistent  with  the  earlier  work  of 
Uberoi  and  Kovasznay  [12].  Corcos  was  then  able  to  show  how  sensor  size 
influenced  resolution  for  high  frequency  fluctuations.  However,  he  was 
not  able  to  give  a  quantitative  estimate  for  the  error. 

Willmarth  and  Roos  [13]  used  an  identical  procedure  with  Corcos  [4] 
to  show  that  Corcos' s  similarity  rule  was  not  valid  for  the  high  fre¬ 


quency  range  when  the  spatial  separation  between  transducers  is  less 
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than  0.78  5*.  Willmarth  and  Wooldridge  [14],  in  their  measurements, 

* 

.  .  .  oj6 

discarded  fluctuations  with  frequencies  less  than  — —  ■  0.14,  because 

Uoo 

of  poor  reproducibility  of  the  power  spectrixn  in  that  range.  They  have 
attributed  the  problem  to  wind  tunnel  noise,  and  propagation  of  disturb¬ 
ances  upstream.  Bull's  [15]  power  spectral  density  has  the  same  drop 

off  and  peak  frequency  as  Willmarth  and  Wooldridge  [14]  but  he  was  able 

.  a)6  * 

to  make  reproducible  measurements  at  frequencies  as  low  as  — —  *  0.025. 


Blake  [6]  used  a  very  small  "pinhole"  microphone  to  study  the  con¬ 
tribution  of  small  scale  fluctuations  to  the  power- spec trura .  He  pointed 

* 

out  that  at  high  frequencies  the  scaling  on  outer  variables  (6  and  U^) 
was  poor.  As  has  been  discussed  earlier,  Bull  and  Thomas  [10]  made 
spectral  measurements  with  two  different  types  of  transducers.  They 
found  that  the  high  frequency  energy  content  was  four  times  greater  when 
measured  with  pinhole  microphones  than  when  measured  with  piezoelectric 
transducers.  However,  by  employing  a  specially  fitted  cap,  the  piezo¬ 
electric  transducer  was  converted  into  a  pinhole  microphone  and  the 
measured  power-spectrum  was  then  in  close  agreement  with  the  spectrum 
produced  by  the  original  pinhole  microphone. 

Bull  and  Langeheineken  [11]  have  pointed  out  the  lack  of  experi¬ 
mental  pressure  data  in  low  Reynolds  nunber  turbulent  flows.  Their  work 
covered  for  the  first  time  a  wide  range  of  Reynolds  numbers  in  a  single 
experiment.  They  have  shown  indeed  that  the  power- spec trura- -particular¬ 
ly  at  high  frequencies — is  highly  Reynolds  number  dependent.  Also, 
their  non-dimensional izat ion  of  power-spectra  (they  have  used  a  variety 
of  dimensionless  correlations)  suggest  that  there  is  no  single  similar- 
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ity  variable  which  causes  the  spectral  curves  to  collapse  onto  a  single 


curve . 


Calculation  of  auto- correlation,  R  (t),  is  very  useful  in  deter- 

PP 

mination  of  the  time  scale  over  which  the  pressure  fluctuations 

correlate  (along  with  determination  of  convection  speed  and  decay  of  the 

pressure  field).  That  is,  the  time  scale  is  indicative  of  the  time 

interval  over  which  fluctuations  "remember"  their  past;  and  the  pressure 

correlation  is  approximately  the  same  as  the  velocity  correlations, 

because  they  are  determined  by  the  same  physical  processes. 

Theoretically  (statistically  is  a  better  word)  it  has  been  understood 

that  any  peak  in  the  power  spectrum  ( <f> (uo ) )  away  from  the  origin  would 

result  in  negative  region  of  R  (t).  However,  the  reverse  does  not  hold 

PP 

true. 


Experimentally  (Bull  [151  and  Blake  [6]),  it  has  been  proven  that 
there  is  a  peak  away  from  the  origin  in  the  power-spectral  density; 
therefore,  one  must  expect  certain  negative  regions  for  the  auto-corre¬ 
lation.  Based  on  his  integration,  Kraichnan  [8]  was  able  to  predict 
that,  unlike  isotropic  turbulence,  any  turbulent  flow  with  a  slight 
anisotropy  will  have  negative  correlation  regions  for  long  separation 
distances.  The  above  suggestion  was  later  confirmed  by  Bull  [15]  and 
Blake  [61. 

Monte  Carlo  methods  have  been  used  in  this  study  to  construct  the 
synthetic  pressure  distributions.  This  work  is  an  extension  of  earlier 
work  by  Ash  [16]  and  Ash  and  Khorrami  [17]  where  detailed  accounts  of 
the  construction  of  the  computer  programs  can  be  found.  Consequently, 
only  the  major  modifications  to  the  previous  work  have  been  discussed 


The  work  which  follows  discusses  simulation  improvements  related  to 


individual  event  construction,  amplification  factors,  amplitude  decay, 
and  wavelength  decay  functions.  Subsequently,  the  statistical  measures 

of  the  quality  of  the  simulation  are  presented  and  discussed.  Finally, 

/ 

an  interpretation  of  the  transducer  size  effect  on  root-mean- square 
pressure  calculations  is  presented  and  conclusions  have  been  made. 
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2.  ANALYSIS 


2.1  Transducers  with  Finite  Size 
The  original  computer  program  (Ash  [16])  generated  pressure  as  an 
instantaneous  point  function  rather  than  the  average  force  which  is 
measured  by  a  transducer  of  finite  area.  The  distance  between  adjacent 
simulation  points  and  the  time  resolution  could  be  specified  by  the 
user,  which  enabled  arbitrary  spatial  resolvability.  However,  in  order 
to  check  the  simulation  and  produce  results  which  could  be  compared  with 
experimental  data,  it  was  necessary  to  introduce  features  which  allowed 
the  influence  of  surface  area  to  be  studied. 

As  described  in  Chap.  1,  the  finite  size  and  the  type  of  a  pres¬ 
sure  transducer  can  create  several  size  effects  (Corcos  [4]  and  Bull  and 
Thomas  [10]).  The  frequency  resolution  of  any  sensor  is  related  direct¬ 
ly  to  its  characteristic  diameter,  as  well  as  its  construction  and  its 
sampling  time.  As  the  size  of  a  transducer  increases,  its  frequency  re¬ 
sponse  decreases.  Furthermore,  a  large  sensor  causes  a  severe  averaging 
effect  (see  appendix  A)  which  can  eliminate  some  pressure  features.  The 
averaging  effect  results  in  erroneous  measurements,  and  causes  an  inac¬ 
curate  representation  of  the  high  frequency  fluctuations. 

With  the  previous  discussion  in  mind,  an  array  of  storage  locations 
(with  arbitrarily  prescribed  size)  were  substituted  for  the  point 
locations  used  in  previous  simulations.  The  program  has  been  modified 
to  permit  arbitrary  specification  of  "sensor"  locations  and  sizes.  The 
span  of  the  "flat  plate"  from  the  front  of  the  first  sensor  to  the  back 
of  the  last  sensor  is  called  the  model  length.  A  schematic  diagram  of 


the  model  is  shown  in  Fig.  1.  The  development  length  is  the  segment 
ahead  of  the  first  sensor  location  which  does  not  contain  a  complete 
simulation  history. 


The  width  of  the  transducer  is  assumed  to  be  equal  to  its  length 
(square  area).  The  range  of  prescribed  diameters  (here  we  use  the  word 

"diameter"  for  convenience)  is  0.0196  <  d  <1.176.  Figure  2  which 

6 


represents  the  spanwise  or  lateral  correlation  data  measured  by 
Willmarth  and  Wooldridge  [14],  indicates  that  the  lateral  correlation 
coefficient  remains  near  unity  over  that  diameter  range.  It  has  been 
assumed  subsequently  that  the  pressure  fluctuations  at  the  center  of  the 
sensing  element  are  not  affected  appreciably  by  adjacent  uncorrelated 
spanwise  fluctuations. 


2.2  Amplification  Factor 

When  a  constant  convection  speed  of  0.8  U—  was  used  in  the  program, 
it  was  found  that  convection  speed  calculated  from  the  space-time  corre¬ 
lation  curves  exceeded  the  free  stream  velocity;  this  is  obviously  in¬ 
consistent  with  every  experimental  finding.  Several  different  correc¬ 
tions  were  tested  before  it  was  concluded  that  the  problem  was  caused  by 
the  instantaneous  creation  of  low  frequency  (large  wavelength)  disturb¬ 
ances.  Some  of  those  individual  disturbance  events  were  created  which 
spanned  more  than  one  simulation  location  at  the  instant  of  creation. 
These  large  eddy  events  would  thus  have  an  infinite  convection  speed. 

Of  course,  the  number  of  such  disturbances  is  small  compared  to  the 
total  but,  since  the  low  frequency  fluctuations  decay  slowly,  the 
convection  speed  effect  can  be  pronounced. 
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An  amplification  time  interval  has  been  employed  that  causes  the 
disturbance  amplitude  to  grow  quickly  from  a  small  "zero"  value  to  full 
amplitude,  in  order  to  eliminate  the  infinite  convection  speed  effect. 


An  equation  of  the  form 


P(x,t)  =  PQ(x,t)  [l 


-  4 ■ 605t2 

(6/U  >2 
c 


has  been  used  where  PQ(x,t)  is  the  randomly  generated  event  amplitude, 

* 

is  the  disturbance  convection  speed  and  6  is  the  displacement 
thickness,  while  t  is  the  relative  time  (0  at  the  instant  of  gener¬ 
ation).  Since  convection  speed  is  frequency  dependent  the  above  growth 
time  interval  becomes  a  larger  fraction  of  the  disturbance  lifetime  for 
high  frequency  fluctuations.  In  fact,  at  very  high  frequencies  the 
amplification  time  can  be  larger  than  the  decay  time,  preventing  the 
fluctuation  from  becoming  significant. 

Figure  3  shows  the  plot  of  the  amplification  function  for  a  fre- 
★ 

quency  of  — =  0.205  which  is  near  the  nominal  peak  frequency.  In 
00 

earlier  simulation  programs,  the  amplification  function  had  been  non-di- 

mensionalized  with  respect  to  boundary  layer  thickness,  5.  However, 

further  testing  has  shown  that  6  was  not  a  proper  choice.  Rather,  the 

★ 

displacement  thickness,  S  ,  scales  the  simulation  results  adequately  to 
enable  close  approximation  of  a  variety  of  experimental  results.  Al¬ 
though  the  amplification  procedure  caused  the  simulated  convection  ve¬ 
locity  to  decrease  below  the  free  stream  velocity,  the  spanning  effect 
due  to  low  frequency  disturbances  still  remained.  That  is,  for  short 


f  .  •  *'  «•  **  mJt  i vV-  -J*.  •-  *.  •.  -k  *, 


separation  distances  the  convection  speed  is  of  order  0.9  UB  compared  to 

an  experimental  value  which  is  of  order  0.6  U  (even  though  the  assumed 

00 

event  convection  speed  is  on  the  order  of  0.6  U  ). 

00 


2.3  Convection  Speed 

Although  the  constant  convection  speed  of  0.8  gave  reasonable 
results,  experiments  have  shown  that  wall  pressure  fluctuations  travel 
downstream  with  frequency  dependent  convection  speed  (5,  15).  In  order 
to  represent  the  experimental  results  more  realistically  it  was 
necessary  to  introduce  a  variable  convection  speed  into  the  simulation. 

A  problem  still  remained,  because  virtually  every  experimenter  has  found 
a  different  range  of  convection  velocities  (see  Emmerling,  Meier  and 
Dinkelacker  [5],  Bull  [15]).  This  discrepancy  between  individual  event 
convection  speed  and  apparent  convection  speed  was  most  pronounced  for 
the  high  frequency  fluctuations  which  are  normally  attributed  to  near 
wall  effects. 

Initially  the  semi-theoretical  convection  speed  model  of  Corcos  [4] 
was  used  to  represent  the  convection  velocity.  His  model  assuned  an 
exponential  variation  in  convection  speed.  However,  in  order  to  avoid 
computational  complexity,  a  linear  variation  of  convection  speed  with 
frequency  was  preferred.  Since  experimental  data  were  so  scattered,  the 
linear  assumption  was  a  reasonable  compromise  between  the  recent  experi¬ 
ments  and  a  simple  functional  form.  The  linearized  convection  speed  was 
assumed  to  decrease  linearly  from  78  percent  of  the  free  stream  velocity 
at  "zero"  frequency  to  a  constant  46  percent  of  the  free  stream  velocity 
at  frequencies  greater  than  or  equal  to  34  times  peak  frequency  in  the 

power  spectrum.  The  peak  frequency  has  been  taken  here  as  f  ■  0.33 

P  1* 
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Hz  from  Bull  [15].  The  assumed  variation  of  convection  speed  with 
frequency  is  shown  in  Fig.  4  along  with  the  limited  experimental  data 
from  Willmarth  and  Wooldridge  [14]. 

2.4  Frequency  Cutoff  Range 

The  upper  limit  of  the  frequency  range  for  individual  events  is 

determined  by  economics  (in  terms  of  computer  time,  and  memory)  and  the 

resolving  time  step,  which  in  turn  is  related  directly  to  the  peak 

frequency.  However,  the  lower  limit  of  the  frequency  was  determined 

from  totally  different  considerations.  First,  in  every  experiment  the 

lower  frequency  fluctuations  have  been  disregarded  because  of  external 

uncorrelated  noise  generation  (due  to  vibration  of  the  wind  tunnel  or 

the  fan).  Since  the  computer  does  not  have  this  problem,  the  low 

frequency  limit  could  be  adjusted  in  agreement  with  flow  physics.  The 

lower  limiting  frequencies  which  can  be  employed  in  this  study  create 

disturbances  which  have  wavelengths  on  the  order  of  50  5,  and  imply 

significant  correlations  over  very  long  distances.  This  is  unrealistic 

in  a  boundary  layer  flow.  For  the  particular  case  considered  here,  the 

(uxS*) 

simulation  generates  frequencies  in  the  Strouhal  number  ~ —  range  of 

0.076  to  15.78,  where  the  lowest  frequency  utilized  corresponds  to  a 
wavelength  of  approximately  8  6 , 

2.5  Amplitude  Decay  Function 

The  earlier  experiments  concerning  the  pressure  fluctuations  have 
revealed  that,  contrary  to  the  Taylor  hypothesis,  eddies  must  decay 
while  they  are  convected  downstream.  Experiments  by  Willmarth  and 
Wooldridge  [14]  and  later  by  Bull  [15]  have  shown  that  each  disturbance 


with  initial  frequency. 


travels  a  distance  which  scales  with  its  wavelength.  An  amplitude  decay 
function  of  the  form: 

AC*.  A)  -A  e-A^A)" 

0 

has  been  employed  where  Aq  is  the  original  amplitude,  X  is  the 
disturbance  wavelength,  x  is  the  distance  from  the  point  of  generation, 
and  a  and  n  are  two  arbitrary  constants,  which  were  determined  by 
trial  and  error.  The  constants  were  specified  ultimately  by  assuming 
that  each  disturbance  travelled  downstream  a  distance  of  five  wave¬ 
lengths  before  it  had  decayed  to  one  percent  of  its  original  amplitude. 

Figure  5  shows  the  variation  of  the  decay  function  with  distance  for  a 
.  .  uiS* 

disturbance  with  a  frequency  of  —  ■  0.205.  The  functional  form  used 

00 

here  is  consistent  with  the  experimental  findings  of  Panton  et  al.  [18] 
which  indicated  the  initial  decay  rate  is  most  rapid.  Also  Fig.  6  shows 
the  combined  effect  of  amplification  factor  and  the  decay  function  on 
the  amplitude. 

2.6  Wavelength  Decay  Function 

The  three  dimensionality  of  a  turbulent  boundary  layer  suggests 
that  the  wavelength  of  a  disturbance  may  change  along  with  its  ampli¬ 
tude.  This  is  a  very  important  concept  when  it  is  thought  of  in  terras 
of  the  power-spectrum  and  auto  correlation.  As  each  disturbance  decays 
in  wavelength,  the  frequency  of  the  disturbance  will  increase  which  in 
turn  shifts  the  peak  frequency,  along  with  the  power  spectrum  curve, 
toward  higher  frequencies. 
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Figure  5.  Assumed  variation  of  amplitude  with  distance  for  a  disturbance  with  an  initial 
frequency  of  ai6*  =  0.205. 


A  wavelength  decay  function  of  the  form 

A(x)  =  A  [k  +  D  e"b(x)  ] 

has  been  used  where  A^  is  the  original,  randomly  generated  wavelength, 
K,  D,  b  and  m  are  arbitrary  constants  which  have  been  adjusted  by 
trial  and  error.  The  wavelength  decay  function  allowed  the  wavelength 
to  be  reduced  to  a  specified  fraction  of  its  original  wavelength  at  the 
time  of  total  decay  of  the  disturbance.  The  wavelength  function  has 
improved  the  simulated  auto-correlation  but  it  cannot  be  considered  an 
optimum  function.  Figure  7  shows  the  variation  of  the  wavelength  decay 
function  with  distance. 

In  a  few  cases  a  frequency  dependent  decay  function  similar  to  the 
one  used  for  amplitude  (see  Sec.  2.5)  was  employed.  However,  that 
frequency  dependent  functional  form  did  not  provide  any  major  improve¬ 
ment  in  simulated  statistics  but  did  create  some  programming  difficulty. 
Therefore,  all  the  runs  presented  in  this  work  have  used  a  wavelength 
decay  of  the  form  given  above. 

To  make  sure  that  a  realistic  signal  was  generated  by  the  simula¬ 
tion,  short  records  of  the  pressure  histories  at  three  different  spatial 
locations  have  been  produced.  Those  histories  are  shown  in  Figs.  8 
through  10,  and  it  is  obvious  that  there  are  recognizable  events  which 
have  been  convected  over  lortg  distances.  However,  there  are  also  sub¬ 
stantial  features  which  represent  "uncorrelated"  information.  There¬ 
fore,  it  was  concluded  that  the  simulated  pressure  records  were  similar 
to  actual  experimental  records. 

The  program  incorporating  the  modifications  described  herein  is 


listed  in  Appendix  B 


Assumed  variation  of  wavelength  decay  function  with  distance 


Simulated  pressure  fluct 


3.  RESULTS  AND  DISCUSSION 


3.1  Two  Point  Correlation 

Two  point  correlations  have  been  calculated  for  separation  dis¬ 
tances,  x  ,  between  1.66  and  13.37  for  a  wide  range  of  transducer 
* 

6 

diameters.  The  auto-correlation  and  bounding  cases  (1.66  and  13.37)  are 
shown  in  Figs.  11  to  13  along  with  the  experimental  data  of  Bull  [15]. 

X 

Figure  11  shows  the  auto-correlation  at  —  =  0. 

6 

It  should  be  noted  that  the  widths  of  the  positive  correlation 

curves  become  narrower  as  the  transducer  diameter  is  decreased.  Note 

also,  that  the  negative  portion  of  auto-correlation  decreases  as  the 

diameter  is  reduced.  At  the  point  where  — —  =  0.0196,  there  is  an 

6 

inflection  in  the  auto-correlation  (at  a  dimensionless  time  delay 

T*0.3).  The  inflection  point  has  been  caused  by  the  resolution  time 

step.  That  is,  at  the  point  where  — —  is  0.0196,  the  time  required  for 

6 

many  of  the  simulated  pressure  events  to  cross  the  transducer  is  smaller 
than  the  user  specified  time  step.  Hence,  the  region  of  large  positive 
correlation  is  diminished  due  to  the  absence  of  those  events. 

A  phenomenon  which  is  similar  to  the  inflection  in  the  auto-corre¬ 
lation  curve  can  be  observed  in  the  two-point  correlations  for  the  small 

separation  distance  case  [— r-  *  1.66,  Fig.  12]  and  for  the  large  separa- 

<5 

tion  distance  case  [Fig.  13].  As  the  transducer  size  is  reduced,  it  is 
capable  of  resolving  smaller  wavelength  disturbances,  but  those 
disturbances  decay  so  rapidly  that  the”  decorrelate  adjacent 


Solid  line  Bull  (15) 

Broken  lines  present  simulation 


Figure  11.  Comparison  between  the  auto-correlation  function  resulting  from  a  simulated  wall 
pressure  history  and  the  experimental  measurements  of  Bull  [151  for  different 


c 

_ ,  0) 

£  $ 
__  Q. 
D  U» 
CD  0* 
C 


2  = 
if  S 

$  I  ■°r« 


S  !S  2 

®  cti  n  »- 


[<•] 


rH 

rH 

03 

* 

u 

T3 

o 

0) 

VM 

4-J 

03 

sO 

rH 

vO 

a 

• 

i 

rH 

•H 

CO 

II 

cd 

x| 

* 

<1 

K? 

e 

0 

u 

j- 

o 

4h 

Mh 

00 

r— i 

c 

m 

•H 

fH 

4-i 

k— 

rH 

3 

rH 

cn 

rH 

a) 

3 

u 

« 

c 

Uh 

o 

O 

•H 

AJ 

CO 

u 

AJ 

c 

c 

3 

0) 

4h 

E 

0) 

c 

V- 

o 

3 

•H 

71 

4H 

cd 

CT3 

0) 

rH 

6 

<U 

rH 

l- 

cd 

0 

•u 

a 

C 

0) 

•u 

e 

C 

•H 

»H 

o 

QJ 

a 

a. 

i 

X 

• 

o 

a) 

CO 

5 

0) 

4J 

a> 

N 

j= 

•H 

0) 

AJ 

CO 

JC 

4J 

T3 

M 

C 

0J 

c 

cd 

a 

0) 

3 

0) 

>N 

"O 

3 

J- 

CO 

aj 

0 

c 

(1) 

aj 

cd 

-Q 

CO 

)H 

fH 

AJ 

C 

x: 

o 

AJ 

<u 

c 

V- 

a) 

»- 

3 

>- 

cd 

CO 

<u 

a 

CO 

Uh 

1 

0) 

4h 

o 

»-< 

•H 

u 

Cu 

*3 

# 

CM 

*— < 

<D 

u 

3 

00 

fH 

U. 

points.  Consequently,  two-point  correlations  are  reduced  as  d_  is 

★ 

6 

decreased  below  0.117. 

A  more  subtle  effect  can  also  be  observed  in  the  two-point  correla¬ 
tion  data  which  is  due  to  improved  resolution  of  the  low  frequency,  long 
wavelength  events.  For  large  transducer  sizes,  a  filtering  effect  is 
created  in  the  opposite  sense  to  the  high  frequency  case  just  discussed. 
The  larger  transducers  cannot  resolve  and,  therefore,  are  not  influenced 
by  the  high  frequency  disturbances.  Consequently,  the  low  frequency 
information  is  less  contaminated  and  produces  a  larger  correlation  coef¬ 
ficient  than  is  actually  present.  This  increased  correlation  effect  is 
manifested  both  in  an  increase  in  maximum  correlation  value  and  an  in¬ 
crease  in  correlation  width  at  a  given  separation  distance. 

Finally,  for  sufficiently  large  transducer  sizes,  the  measured 
pressure  signal  becomes  contaminated  by  the  splitting  effect  discussed 
in  Appendix  A.  That  is,  a  continuous  pressure  event  is  split  into  two 
events  separated  by  a  "dead  band"  when  the  particular  pressure  event  re¬ 
sides  entirely  upon  the  transducer.  This  splitting  effect  causes  a  de- 

correlation  to  occur  between  adjacent  points  which  evidently  begins  for 

d 

transducer  sizes  on  the  order  of  — r  ■  0.117  in  this  case. 

6 

Filtered  correlation  data  have  been  reported  by  Willmarth  and 

★ 

Wooldridge  (14)  for  pressure  signals  in  the  0.41  <  — —  <  0.95  (low  band) 

Uoo 

★ 

and  4.1  <  — — <  6.8  (high  band)  frequency  ranges.  In  order  to  compare 
00 

the  correlation  data  produced  by  the  simulation  with  these  results,  a 
pseudo  band  pass  filter  has  been  employed.  A  comparison  can  be  made  by 


employing  IF  STATEMENTS  in  the  simulation  program  rather  than  digital 
filters.  That  is,  only  those  randomly  generated  pressure  events  which 
had  frequencies  in  the  prescribed  intervals  were  allowed  to  contribute 
to  the  transducer  histories.  Since  there  are  higher  harmonics 
associated  with  the  discrete  pressure  events,  (which  were  allowed  to 
accumulate  using  this  approach),  the  "filtering"  was  not  perfect. 
However,  the  actual  experiments  did  not  employ  perfect  filters  either. 
The  low  band  and  high  band  data  are  shown  in  Figs.  14  and  15  along  with 
the  experimental  data  of  Willmarth  and  Wooldridge  (14).  In  Fig.  15, 
only  the  locus  curve  of  the  space-time  correlation  peaks  is  shown  as  a 
solid  line.  The  apparent  convection  speeds  for  the  simulated  case  are 
approximately  0.6  U  for  the  high  band  and  0.8  for  the  low  band. 

Narrow  band  correlation  data  have  also  been  produced  using  a 


similar  pseudo  filter  for  a  spatial  separation  of 


5.0.  The  fre¬ 


quency  bands  were  defined  as  16.0  percent  on  either  side  of  centered 


frequencies  of  -jr—  =  0.33,  0.62,  1.15  and  2.39  and  are  shown  in  Fig. 

16.  The  narrow  band  convection  speeds  show  the  proper  celerity  de¬ 
crease  with  increasing  frequency.  A  center  frequency  of  0.14  was 


also  examined  but  the  correlation  results  were  unsatisfactory  for  rea¬ 
sons  which  will  be  discussed  in  the  next  section. 


3.2  Power  Spectra 

Power  spectra  for  the  different  sized  transducer  simulations  can  be 


developed  directly  from  the  auto-correlation  data.  However,  interpret  a- 


Willmarth  and  Wooldridge  (14) 


Figure  14.  Filtered  longitudinal  space-time  correlations,  low-frequency  band 


<6.8  present  work 


Filtered  longitudinal  space-time  correlations,  high-frequency  band 


tion  of  power  spectral  data  is  limited  to  a  frequency  interval  which  is 
controlled  by  the  sampling  data  (specified  time  step)  and  the  total  time 
record.  Since  either  a  Fourier  transform  or  a  fast  Fourier  transform 
must  be  used  to  extract  the  power  spectrum  from  the  auto-correlation, 
consideration  must  be  given  to  the  upper  frequency  limit  (Nyquist  fre¬ 
quency)  where  the  inverted  Fourier  transform  will  begin  repeating 
itself.  The  sampling  frequency  f  must  be  greater  than  or  equal  to  2  f^ 
(where  f  is  the  Nyquist  frequency)  in  order  to  be  able  to  recover  the 
original  signal  from  the  Fourier  transformation. 

Figures  17  through  20  represent  the  variation  of  povrer  spectra  with 
transducer  diameter,  where  the  shaded  area  is  the  range  of  experimental 
data  given  by  Bull  [15],  Willmarth  and  Wooldridge  [14],  and  Blake  [6], 
The  maximum  and  minimum  frequencies  generated  in  the  simulation  are 

f  <  1 

max  7— 

At 

where  A t  is  the  resolving  time  step  and 

f  .  >  1 

min  — 

where  T  is  the  total  record  time.  The  above  maximum  frequency  gener¬ 
ated  by  the  simulation  is  smaller  than  the  sampling  frequency  f  associ- 

s 

ated  with  Nyquist  frequency  f^  of  turbulent  wall  pressure  fluctuations. 
Therefore,  the  discrepancies  between  the  experimental  data  and  the  simu¬ 
lated  spectra  at  the  high-  and  low-frequency  extremes  are  due  to  the 
time  step  resolution  and  the  shortness  of  the  simulated  record,  respect¬ 
ively.  However,  the  results  are  in  close  agreement  with  the 
experimental  data  in  the  raid- frequency  range  for  a  wide  variety  of 
diameters.  The  above  figures  indicate  a  small  increase  in  the  power 


level  at  high  frequencies  as  the  diameter  approaches  zero.  Considering 


the  degree  of  approximation  used  in  the  present  work  and  the  scattering 
of  experimental  data  the  agreement  between  the  simulation  and 
experimental  results  is  considered  satisfactory. 


3.3  Variation  of  P 

R.M.S. 

For  many  years  experimenters  scaled  transducer  dianeter  data  with 

★ 

the  boundary  layer  displacement  thickness  6  [13-15]  and  showed  only 

small  variations  of  P  as  the  parameter  d  approached  zero.  How- 

R.m.s.  — 

6 

ever,  more  recent  results  suggest  that  small  scale  pressure  fluctuations 
in  turbulent  boundary  layers  may  have  a  significant  effect  on  the 
P  „  value  (Blake  [6],  Emerling,  Meier  and  Dinkelacker  [5],  and  Bull 
and  Ihomas  [10]).  Since  most  of  these  high  frequency  fluctuations  are 
produced  in  the  wall  region,  it  was  suggested  by  Emmerling,  Meier  and 
Dinkelacker  [5]  that  the  scaling  of  transducer  diameter  should  be  with 
respect  to  the  inner  variable  _v  which  expands  the  lower  limiting 


region  significantly.  They  also  suggested  that  instead  of  wall  shear 
stress,  free  stream  dynamic  pressure  should  be  used  in  normalizing  the 
R.M.S.  pressure. 

For  the  present  study,  it  was  decided  to  concentrate  numerical 

testing  in  the  range  <  200,  in  order  to  verify  the  above  sugges- 

v 

U  d 

tion.  Therefore,  diameters  with  the  value  — —  -  9,  18,  35,  54,  72,  90, 
135,  180,  271,  360,  451,  and  541  have  been  used  for  testing.  The  ampli¬ 
tude  of  P  can  be  adjusted  in  this  simulation  by  simply  changing  a 

scaling  parameter  in  the  program.  However,  an  adjustment  was  made  only 
for  one  large  diameter  case  to  establish  a  realistic  level  for  use  with 


all  other  cases.  Hence,  each  subsequent  case  was  subjected  to  an 
identical  simulated  pressure  history.  Figure  21  shows  the  variation  of 
Pr  m  s  with  the  transducer  size  as  measured  by  Emmerling,  Meier  and 
Dinkelacker  [5]  and  Bull  and  Thomas  [10]  along  with  the  present  simula¬ 
tion  results.  Although  there  is  a  fluctuation  of  P  ,  values  for 

R .  M .  S . 

U  d.  . 

t  >  200,  the  oscillation  is  small  and  P  can  be  represented  by  a 

K«  Mi  b  • 

straight  line  for  large  diameters.  The  constant  simulated  P 

R . M. S . 

region  is  also  followed  by  a  sharp  increase  in  the  R.M.S.  value  as  the 
diameter  shrinks.  The  simulated  variation  approximates  the  curve  given 
by  Bull  and  Thomas  [10]. 

Although  the  increase  in  P  is  significant,  it  was  not  clear 

whether  the  jump  was  caused  by  the  inclusion  of  high  frequency 
disturbances  or  by  achieving  a  better  resolvability  of  large  disturb¬ 
ances  at  small  diameters.  Therefore,  a  simulation  was  conducted  in 

V 

which  all  the  disturbance  events  with  wavelengths  smaller  than  — -  *  180 

were  excluded  from  the  simulation.  The  remaining  fluctuations  were 

U  d 

resolved  with  a  diameter  of  — —  =  18  which  is  in  the  range  of  increased 
P„  „  _  .  That  test  produced  the  same  P  and  has  indicated  that 

very  high  frequency  pressure  fluctuations  do  not  contribute  measurably 


to  the  R.M.S.  level. 


4.  CONCLUSION 


A  wide  range  of  transducer  diameters  have  been  used  in  the  present 
simulation  to  study  the  variation  of  R.M.S.  pressure  with  sensor  size 
when  the  transducers  are  exposed  to  identical  signals.  It  has  been 
shown  that  the  P  c  level  increases  moderately  at  small  values  of 

t\i  Hi  b  i 

U  d 

- —  which  is  consistent  with  more  recent  experimental  data.  The 

simulated  result  indicates  that  P  „  increases  from  a  nominally 

R. M.  • 

UTd 

constant  value  for  sensor  sizes  smaller  than  -  =  40  which  differs  from 

v 

U  d 

the  departure  point  of  — —  =  170  suggested  by  the  experimental  data. 

The  simulated  curve  closely  approximates  the  result  given  by  Bull  and 
Thomas  [10]  in  terms  of  the  magnitude  of  increase. 

Although  the  rise  in  the  R.M.S.  level  is  significant,  the  result¬ 
ing  curve  suggests  that  the  increase  is  not  as  large  as  that  reported 
earlier  by  Emraerling,  Meier  and  Dinkelacker  [5].  It  is  known  that  the 
above  contribution  to  Pr  ^  g  has  been  attributed  to  small  scale 
structures  in  the  inner  region  of  the  turbulent  boundary  layer  by 
several  workers  in  the  field.  However,  the  R.M.S.  level  of  the  simulat¬ 
ed  signal  which  deleted  high  frequency  events  indicated  that  the  contri¬ 
bution  of  small  scale  fluctuations  to  the  P_  _  is  insignificant. 

Power-spectra  have  been  calculated  for  several  different  transducer 
diameters.  The  power  spectra  showed  slight  increases  in  energy  content 
at  high  frequencies  and  a  slight  increase  in  the  peak  frequency  as  the 
diameter  was  reduced.  The  shapes  of  the  spectra  are  very  similar,  and 
the  agreement  with  experimental  data  is  reasonably  good. 
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Two  point  correlations  have  been  calculated  for  separation  dis- 

x 

tances  of  —  =  1.66,  and  13.37.  The  correlations  show  clearly  the 

6 

convection  and  decay  patterns  of  the  pressure  field.  For  large  trans¬ 
ducer  diameters  the  width  of  the  correlation  curves  increases.  This  is 
in  accordance  with  the  expectation  that  higher  correlation  is  obtained 
when  the  resolvability  of  small  scale  structures  is  poor.  As  the 
diameter  is  reduced  the  correlation  curves  become  narrower.  Again  this 
is  expected  since,  the  inclusion  of  small  pressure  fluctuations  which  do 
not  correlate  over  large  separation  distances  or  time  delays,  have  a 
significant  decorrelation  effect.  Low-  and  high-frequency  band  correla¬ 
tion  of  the  generated  pressure  fluctuations  also  showed  the  correct 
convection  speeds.  The  narrow  band,  two  point  correlations  showed  a 
decreasing  convection  speed  as  the  center  frequency  was  increased. 
Overall  the  agreement  between  the  simulated  curves  and  the  experimental 
data  of  Bull  [15]  and  Willmarth  and  Wooldridge  [14]  is  quite  good. 

In  summary,  we  have  shown  that  a  simulated  turbulent  wall  pressure 
signal  with  properties  similar  to  experimental  measurements  produces 
root  mean  square  pressure  levels  which  vary  with  sensor  size  in  a  manner 
which  is  similar  to  experiment.  The  size  effect  is  not  pronounced  for 
the  power  spectra,  but  the  interesting  result  that  the  maximum  values  of 
two-point  correlations  decrease  with  decreasing  size  was  observed.  This 
effect  can  be  attributed  to  the  increasing  contribution  from  small 
wavelength  pressure  fluctuations. 
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APPENDIX  A 


FORCE  CONTRIBUTION  OF  A  SINGLE  PRESSURE  EVENT 


The  streamwise  length  of  the  transducer  (size)  and  the  sampling 
rate  (time  step)  both  distort  the  recorded  pressure  signal.  Depending 
on  the  wavelength  (and  frequency)  of  the  particular  fluctuation  event, 
the  net  pressure  contribution  sensed  by  a  transducer  can  vanish  even 
though  the  disturbance  is  over  the  transducer.  Also,  if  sampling  rates 
are  slow,  an  event  may  be  convected  completely  across  the  transducer  in 
less  than  one  time  step.  In  order  to  understand  the  distortion  produced 
by  finite  pressure  transducers  and  sampling  rates,  calculations  have 
been  made  for  a  variety  of  prototype  pressure  events,  AP(x,t),  given  by: 


AP(x,t)  * 


0 

sin  2^(x-ct)  -  1/2  sin  4ir(x-ct)  ,  ct-L 
L  L 


x  <■  ct-L 
<  x  c  ct 


0  ,  x  >  ct 

(t  >  0) 


where  L  is  the  wavelength  and  c  is  the  convection  speed.  Each  event  was 
assumed  to  have  unit  magnitude  and  no  decay  of  any  type  was  assumed  in 
order  to  show  clearly  the  spatial  and  temporal  resolution  effects. 

Since  the  wave  front  for  AP(x,t),  designated  Xf(t),  is  located  by: 


a  contribution  to  the  pressure  force  on  a  transducer  located  on  the  in¬ 
terval,  Xj  <  x<  X2 ,  will  be  made  when 

Xf(t)  >  xx 

and  X  t)  _  L  <  X2  . 

where  the  origin  of  the  coordinate  is  sitting  on  the  wavefront  and  is 

moving  with  the  event  convection  speed,  c. 

Assuming  the  spanwise  (Z)  dimension  of  the  transducer  is  unity,  the 

incremental  force  contribution  due  to  this  pressure  event  at  some  time, 

tn,  which  satisfies  the  wave  front  criteria  (x^  <  t  <  X2+L)  is  given 

c  c 

by: 


L  ©2 

F(tn) .=  ~  J"  [sin  9  -  1/2  sin  20]d9 


cos  29o  -  cos  20 


2"F 


^  +  cos  0i  -  cos  02  =  - £" 


where 


2ttxi  when  Xg  -  L  <  xi 


2*  (X  -L)  _ 

_ f  when  X^  -  L  >  Xi 


and 


02 


2TTX 


f  ,  when  X^  <  x2 


2ttxq  ,  when 
L 


Xf  >  x2 . 


49 


Consequently,  the  average  incremental  pressure  contribution,  AP^,  can  be 
written: 


n  x2_xi  2*  (x2  "  xi 


y  1 1/4  [cos  202_cos  20^]  +  cos  0^-  cos  0, 


The  time  step  discussion  will  be  deferred  temporarily  and  attention 

will  be  restricted  to  some  time,  t  ,  when  the  disturbance  is  over  the 

n 

transducer.  For  very  early  times,  the  wave  front  will  be  over  the 
transducer,  while  the  trailing  edge  of  the  wave  has  not  yet  crossed  the 
boundary,  xi  .  (See  Fig.  Al).  Then, 


while, 


2ttx 

•i-r’ 

so  that 


n  2rr  (x^x  ^ 


1/4  [ 


4ttx^  2ttx^  2ttc1 

-  cos  — —  ]  +  cos  — - cos  — — 


For  very  late  times, 


2ir(ct  _L) 


while, 


♦Integration  can  be  checked  by  taking  the  limit  as  x1  -*■  X2>  which 
restricts  time  so  that  0?  0i  . 


prototype  double  sine  wave  moving  over  a  transducer 


leaving 


ao  -  t  i  /,  r  4itx.  2irct  I  ,  2itct  2itx„ 

iPn  =  L  1/4  [  cos  _ 2_  -  cos  _ nj  +  cos  _ n  -  cos  _ 2_ 

U  2jt(x2~x^)  L  L  L  L 

Two  types  of  wavelength  cases  are  delineated  by  the  case  when 
=  L.  When  the  wavelength  is  greater  than  the  transducer  length, 

L  >  x^-x^  ,  there  will  be  times,  tR,  when  neither  the  wave  front  nor  its 
trailing  edge  will  be  over  the  transducer,  in  which  case 


The  second  case  in  which  L  <  x^-x^,  is  discussed  next. 

Figure  A 2  represents  the  averaged  force  as  a  function  of  time  for  a 
large  transducer  (L  <  x^'X^).  That  figure  shows  how  the  shape  of  the 
double  sine  wave  is  reinterpreted  such  that  the  negative  and  the  posi¬ 
tive  pressure  regions  are  separated  by  a  dead  band.  This  dead  band  is 
created  during  the  time  interval  when  the  entire  wave  is  over  the  trans- 
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ducer.  The  same  effect  is  also  responsible  for  a  significant  change  in 
the  wavelength  of  the  averaged  transient  force.  The  indicated  wave¬ 
length  will  be  twice  that  of  the  event  when  the  transducer  length, 

X2~X1*  ecJual  to  **•  When  greater  than  L,  the  apparent  wave¬ 

length  is  increased  by  more  than  a  factor  of  two,  as  shown  in  the  fig¬ 
ure.  Also  note  that  the  averaging  caused  by  the  transducer  has  a  severe 
effect  on  the  amplitude  of  the  wave. 

At  this  point,  it  should  be  pointed  out  that,  in  order  to  show  the 
above  effects  as  clearly  as  possible,  the  smooth  curve  of  Fig.  A2  has 
been  obtained  by  employing  a  resolving  time  step  which  was  an  order  of 
magnitude  smaller  than  the  time  step  used  in  the  simulation.  Therefore, 
the  same  case  was  repeated  employing  the  nominal  simulation  time 
increment  used  in  this  study.  The  resulting  curve  is  shown  in  Fig.  A3. 
Although  it  is  not  a  smooth  curve  anymore,  the  major  features  described 
above  still  are  detectable. 

A  series  of  tests  have  been  run  for  the  case  in  which  L  >  x^-x^. 

In  those  tests  the  width  of  the  sensing  area  has  been  kept  constant  at  a 
d 

value  of  — *  =  0.176,  while  the  frequency  of  the  wave  and  the  resolving 
6 

time  step  have  been  varied.  Figure  A4  shows  a  typical  wave  with  a 
frequency  of  300  Hz  (Sfc  *  0.287)  resolved  by  a  time  step  which  is  an 
order  of  magnitude  larger  than  the  time  step  used  in  the  simulation. 

The  solid  lines  are  the  numerically  integrated  values,  while  the  squares 
represent  the  theoretical  wave  at  a  point  location.  Although,  the 
agreement  between  the  two  types  of  calculations  is  good,  the  figure  does 
not  carry  sufficient  information  to  represent  the  prototype  double  sine 
wave.  Figure  A5  shows  the  same  wave  resolved  by  a  time  step  equal  to 
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Averaged  force  as  a  function  of  time  with  a  time  step  equal  to  the  simulation  time 
step  (At) . 


□  THEORETICAL  VALUE 
—  INTEGRATED  VALUE 


HI 

3  LU 
_J  Z> 

<  -I 

** 

<  Q 

CJ  LU 


Temporal  resolution  of  a  large  wavelength  double  sine  wave  with  a  time  step  equa 
to  the  simulation  time  step  (At). 


the  simulation  time  step.  The  agreement  between  the  two  methods  is  ex¬ 
cellent  and  the  shape  of  a  double  sine  wave  is  retained.  Subsequently, 


the  frequency  of  the  wave  was  increased  to  2500  Hz  (St=2.39)  while  the 
resolving  time  step  was  kept  constant.  Figure  A6  suggests  that  as  the 
frequency  of  the  wave  is  increased,  the  temporal  resolution  becomes  a 
problem.  Also,  note  that  there  is  a  phase  shift  between  the  two  types 
of  calculations.  Figure  A7  represents  a  wave  with  a  frequency  of  2500 
Hz  resolved  by  a  time  step  which  is  an  order  of  magnitude  smaller  than 
the  time  step  used  in  the  simulation.  The  above  figure  indicates  that 
the  resolution  problem  disappears  as  the  time  step  approaches  zero. 

In  summary,  significant  alterations  of  prototype  high  frequency 
pressure  fluctuation  events  can  occur  due  to  transducer  size  and  resolv¬ 
ing  time  increment.  But,  as  the  wavelength  of  the  individual  event  in¬ 
creases,  the  effects  of  transducer  size  and  temporal  resolution 


diminishes. 


□  THEORETICAL  VALUE 
—  INTEGRATED  VALUE 


APPENDIX  B 


COMPUTER  PROGRAM  LISTING 


PROGRAM  SIMU< INPUT* OUTPUT >TAPF5=  INPUT, TAP E6* OUT  PUT* TAPE 1*  TAPS 2 
1*  TAPE7*T APF8 ) 

DIMENSION  P  (  5* 16500) 

DIMENSION  SUMX(8000)*SUMY< BOOO) 

EQUIVALENCE ( P < 60000 ), S UMX ) , ( P < 7 4000 > *  SUMY  ) 

DIMENSION  PA (  6  )  *  PR {  6  )*  STR T ( 6 ) *  ST  P ( 6 ) 

NR  ■  5 
NW-6 

NDIM-MAXIMUM  TIME  DIMENSION 
NDIM-16500 
N2C-NDIM/500 
CDIM-NDIM 

COR X  IS  THE  ADJUSTMENT  FOR  DX  STREAK. 

IF  MORE  THAN  ONE  BURST  PER  STREAK*  CORX  IS  NOT  UNITY 
CORX-1.5 

CORT  IS  THE  TIME  ADJUSTMENT  FOR  DT-STREAK 
CORT-2. 

XD-DEVELOPMENT  LENGTH  (M) 

XM-MODEL  LENGTH  (M) 

PAR  IS  THE  DECAY  ADJUSTMENT  ACCOUNTING  FOR  DISTURBANCE  SIZE 
PAR-10. 

XD-4.145 

D XT  IS  THE  DISTANCE  BETWEEN  STORAGE  LOCATIONS 

DXT-0.008 

DWT-DXT 

DWT  IS  THE  WIDTH  OF  THE  PRESSURE  STORAGE  LOCATION 

DWT-0.002 

IMJ-0 

XM-0.0 

157  XD-XD+XM 
XT-XD+XM 

DXT-NODE  SPACING  ON  MODEL(M) 

NXM  is  THE  NUMBER  OF  STORAGE  LOC AT  IONS . . . YOU  MUST  CHANGE  THE 

NUMBERS  IN  THE  DIMENSION  STATEMENTS  AND  IN  THE  TAPE  OUTPUT 

NXM-3 

STRT ( 1 ) »XD 

STP(1 )  *  STRT (1 )  +  DWT 

STRT(2)-XD+0.00B466 

S  TP  <  2 ) -STRT(2)+DWT 

STRT (3)«XD+0.06B1B7 

STP(3)-STRT(3)+DWT 

XSMX-2.0+STP (NXM) 

XSMX  IS  THE  RIGHT-MOST  LOCATION  STORING  PRESSURE 

US-33.5 

BLT-0.040B 

UFRIC-1.25B7 

CNU-O. 00001395 

RHO-1.  2 

DPT-BLT/B 


o  o  o  o  o 


C  ENTER  OTHER  DISPLACEMENT  THICKNESS  HERE  IF  DESIRED  (CM) 

C  DIMENSIONLESS  STORAGE  LOC.  W IDTH — W/DISP .  THICK. 

WDM«OWT/DPT 

tf R ITE ( NW»  205 )  WDM 

205  FORMAT( ///>5X>*DIM.  STORAGE  WIDTH  EO.  *>E16.6>//) 

WRITE(NW>206) 

206  FORMAT ( 5X»*L0CATI0N  NO.  DIM.  LOC.*>/> 

DO  30  I *  1 >  NX  M 

XDM«$TRT( I J/DPT 
WR  IT  E  (  NW>  210  )  I »  XDM 

210  FORMAT (10X*I4*10X*  E16. 8) 

30  CONTINUE 

C  CALCULATION  OF  WALL  SHEAR  STRESS*  TW 

TW«RHO*UFRIC*UFRIC 

C  CALCULATION  OF  RMS  PRESSURE  FLUCTUATION 

PRMS*2. O+TW 

C  CALCULATION  OF  MINIMUM  RELEVANT  DISTURBANCE  LENGTH 

C  TWV  IS  TIME  DURATION  OF  DISTURBANCE 

TWV»DPT/US 

C  LENGTH  OF  SMALL  DISTURBANCE 

CLM»0.7*US*TWV 

C  DECAY  BIAS  PARAMETER 

CWV«PAR*CLM 
C  WV*  0 • 007 

C  MISCELLANEOUS  CONSTANTS 

PI*3. 1A15926 
TPI-2.*PI 
HPI«PI/2. 

C0«2. 515517 
C1*0. 802P  53 
C2«0. 010328 
Dl»l .A32788 
D2«0. 189269 
D3«0. 001308 

C  STARTER  FOR  RANDOM  NUMBER  GENERATOR 

XSTART«77294. 

RNM»URAN(XSTART) 

XSTART  «0.0 

C  CALCULATION  OF  NOMINAL  PEAK  FREQUENCY 

FPEAK«0.2057A*US/(TPI*DPT> 

FMAX«10.*FPEAK 

C  COMPATABLE  TIME  STEP 

DTT«1./FMAX/10. 

A  USER  SPECIFIED  TIME  STEP  CAN  BE  ENTERED  HERE 

THIS  TIME  STEP  IS  THE  TIME  INCREMENT  USED  IN  THE  RESOLUTION  OF  THE 
OUTPUT— THE  INTERNAL  FLUCTUATION  TIME  STEP  IS  RANDOM 
MAXIMUM  TIME  IS  CONSTRAINED  BY  COMPUTER  STORAGE.  SET  THE  NUMBER 
OF  ALLOWABLE  TIME  STEPS— NTM 
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NTM-16500 

C  INITIAL  TIME  IS  TO(SEC) 

TO-O. 

DTSUM-O. 

NCT-1 

C  PREFIXES  FOP  RANDOM  NUMBER  CALCULATIONS 

PX«DPT*UFRIC/US 
PT-DPT/US 
PW-US/DPT 

C  CALCULATION  OF  REQUIRED  START  UP  TIME  FOR  SIMULATION 

TS0«1.8*XT/US 
NMIN-TSO/DTT 
TMAX1«NMIN*DTT 

C  IF  NM IN  IS  GREATER  THAN  NTM,  NTM  IS  OVERRIDDEN 

IF(NMIN-NTM)201,  201#  202 
202  NTM-NMIN 
201  CONTINUE 

TMAX«NTM*DTT 

WRITE ( 6# 100)  DTT  ,TMAX 

100  FORMAT <5X,*DT-*, FI 0.8, 10X,*TMAX-*,F10.4,//> 
TMAXS«CDIM*DTT 
TMAX-TMAX+TMAX1 
NTM2 “NO IM— NM IN+2 
TREF-TMAX1 
NREL-0 

C  INITIALIZE  PRESSURE  ARRAY 

DO  52  NX-1, NXM 
DO  52  NT-1, NDIM 
52  P(NX,NT)»0.0 
NFLG-0 
TSUB-O. 

C  INITIALIZE  LOCATION  AND  TIME  BASE,  ETC. 

1  X-0. 

DTAVG-DTSUM/NCT 

NCT-0 

DTSUM-O. 

TO-TO+DT  AVG 
IF(T0-TMAX1)2,150, 150 

150  IF(NFLG)  151,151,15 

151  NFLG-1 
TMAXl-TMAX 

DO  153  J«NTM2,NDIM 
P(I,J)*0#0 
153  CONTINUE 
TO-O. 

15*  CONTINUE 
TSUB-TREF 
TMAXS-TMAXS^TREF 
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2  RNM-URAN(XSTART) 

PNM-O. 005+0. 99*RNM 

C  CALCULATION  OF  DX  USING  RANDOM  NUMBER  RNM 

HPII-HPI+RNM 

DX«PX*(32.2-2/(RNM+0.0619)+72.*RNM**2+0.63*TAN(HPII) ) 

DX*CORX*DX 

X-X+DX 

RNM-URAN(XSTART) 

RNM-0.04+0.955*RNM 

C  CALCULATION  OF  RADIAN  FREQ.  FROM  NEW  RNM 

SRNM-RNM**0.6667 
RRNM«(1.-RNM)**0.74 
FRNM-1 / RRNM—1 

W«PW*(0.523*FRNM+0.799*SRNM-0.785*RNM) 

F-W/TPI 

IF(F.GT.7300. )G0  TO  85 


B5 

87 


C 

C 

C 

C 


C 

C 

C 

3 

971 


961 

962 

963 
43 

964 
C 


DC-0.78-F/ (105. *F PEAK) 

UC«DC*US 
GO  TO  87 
UC-0 .46*US 
TP-l/F 

BOBF AC -4.605/ ( (DPT/UC ) **2) 

DXE-UC*TP 

XQ-X-DXE 

XO  IS  THE  ORIGIN  OF  THE  SINE  WAVE  FLUCTUATION 

X  IS  THE  FRONT  OF  THE  SINE  WAVE 

XS  IS  THE  FIRST  STATION  AT  WHICH  P  IS  RECORDED. 

XS  AND  N  X I  WILL  BE  TAKEN  AS  THE  FIRST  STATION  VALUES — THEN  OVER  I  DDE 

XS-XD 
NX  I-  1 

CHECK  TO  SEE  IF  THE  DISTURBANCE  IS  OVER  THE  MODEL 
I F ( X  D-X ) 3#  3#  5 

IF  THE  DISTURBANCE  IS  OVER  THE  MODEL#  HAS  IT  PASSED  THE  LAST  DATA 
STATION 
CONTINUE 

IF(X0-XSMX)971#1#1 

CONTINUE 

IF (XO.LT.O. )G0  TO  964 
DO  43  I* 1 # NX M 
NNP- XO/STP{I  ) 

IF(NNP)961#962#963 
NX  I  -1 
GO  TO  964 
NX  I  - 1 
GO  TO  964 
CONTINUE 
CONTINUE 
CONTINUE 

N X I  IS  THE  NUMBER  OF  THE  FIRST  STORAGE  LOCATION 


CXNI»NXI 

CONTINUE 

GENERATION  OF  RANDOM  TIME  STEP 
RNM«URAN( XSTART) 

RNM«0. 005+0. 99*RNM 
HPII>HPI*RNM 

DT»PT*( 32. 2-2./  ( RNM+O.  0619 )  +72 .  *RNM**2+0 . 63*T  AN  (  HP  1 1 )  ) 

DT«CORT*DT 

T-TO+DT 

NCT«NCT+1 

DTSUM*DTSUM+DT 

IF(XO.GT.STP(NXM))GO  TO  2 

GENERATION  OF  GAUSSIAN  RANDOM  PRESSURE  AMPLITUDE 
RNM«URAN(XSTART) 

ClND-RNM+0.5 
I ND*  C I ND 
CIND«IND 

PPP*2.*(1.-CIND)-1. 

ARGR*RNM/( 1.+CIND) 

ARG-1./ ( ARGR*ARGR) 

CT* ALOG ( ARG ) 

CM*SQRT(CT ) 

PMG«CM-(C0+CM'MC1  +  CM*C2)  )/(l+CM*(Dl  +  CM*(D2  +  CM*D3)  )  ) 

XS-XS-DXT 

PPP-5.9000+PPP 

PE»PRMS*PMG*PPP 

PE»2.0*PE 

MODIFICATION  OF  PRESSURE  AMPLITUDE  TO  ADJUST  RMS  PRESSURE  LEVEL 
DO  LOOP  FOR  STEPPING  THROUGH  MODEL  STORAGE  LOCATIONS 
DXXE*DXE 

DO  1A  NX«NXI>NXM 
DXE«DXXE 

MODEL  STATION  X-LOCATION 
XS*STRT (NX) 

ARRIVAL  TIME  OF  PRESSURE  FLUCTUATION 

DXS*XS-X 

ADX* A8S ( DXS ) 

DXS«(DXS+ADX )/2. 

TGO-DXS/UC+T 

FLUCTUATION  DEPARTURE  TIME 

DXO«STP(NX )-X0 

TSTP-DXO/UC+T 

DOFS  TST P  EXCEED  TMAX 

IF(TSTP-TMAXS )9,9, 7 

IF(TG0-TMAXS)8#8*13 

NSTOP-NDIM 

GO  TO  10 

NSTOP«TSTP/DTT 

NGO-TGO/DTT 


IF(NG0-NREL>99, 99,98 
98  CONTINUE 

AGX»DX0**0.6 
ARG*-4 , 20* AGX 
DH»0.35+0.65*EXP{ ARG) 

OX  EE* OH* OX E 
OR*OXt-DXEE 
R  E*OR/ ( DTT*UC  ) 

RE«RE+0.5 

NP*RE 

NSTOp*NSTOP-NP 

DXE-DXEE 

IF(DXE. IE. 0.0001 )G0  TO  2 

C  00  LOOP  FOR  SUCCESSIVE  TIME  CONTRIBUTIONS  TO  THE  SAME  X  LOCATION 

00  12  NT*  NGO, NSTOP 
TC«NT*OTT 
DELT*TC-T 
XTF*  X  +  UC* (CELT) 

XTB  *  XT  F-OX  E 

THET1*(XTF-STP(NX) )*TPI/DXE 
TH1A-ABS (THET1) 

THET1* ( THET1+TH1 A ) /2 . 

THET2*(STRT(NX)-XTB)/0XE 
TH2A«ABS (THET2 ) 

THET2»TPI-(THET2+TH2A)*PI 

XOT«UC*DELT 

IF (XOT-O. 0000001)19, 19, 18 

18  CONTINUE 

C  FREQUENCY  DEPENDENT  VARIABLE  DECAY  RATE  ADJUSTMENT 

AR-XOT/OXXE 
ARGX«-2.8A*(AR**0.30) 

IF(ARGX.LE.-30. )G0  TO  73 
DECA«EXP ( ARGX  ) 

GO  TO  23 
73  OECA-O.O 
GO  TO  23 

19  CONTINUE 
DEC A«1 . 

23  CONTINUE 

DTH1*  2 • *THET1 
DTH2*  2 . *THET  2 

CF1«(C0S(DTH2)-C0S(DTH1) )/A. 

CF2«C0S (THETl)-COS (THET2) 

DP*  PE  * ( CF 1+CF  2 ) 

DP«DP*DECA*DXE/ (TPI*DWT) 

B0BBA«B0BFAC*DELT**2SIF (B0BBA.GT.300. )BOBBA-300.0 
DP«DP*( l.-EXP(-BOBBA) ) 

NTIME-NT-NREL 

P(NX,NTIME)«P(NX,NTIME)+OP 
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12  CONTINUE 
99  CONTINUE 
13  CONTINUE 
1 A  CONTINUE 

130  FORMAT (5X>F10.6*6E14.6) 

GO  TO  2 
15  CONTINUE 
NURT-NTM/8 
CNT-NTM 
DO  11  J “1  *  NX  M 
P  A  (  J  )«0.0 
11  P  R  (  J  )  »  0 . 0 

DO  21  .  J«1,NXM 
DO  20  I«l* NTM 
PFI«P(J,I) 

PA ( J ) ■ P A  C  J  )  +PF I 
P  R ( J ) ■ PR ( J )  +  PFI*PFI 
20  CONTINUE 

PA(J)«PA(J)/CNT 
WRITE(NW,192)  Jj.PA(  J) 

PSJ-PA ( J ) *PA ( J  ) 

PRJ-PR( J) /CNT 
PR(J)*SQRT (PRJ-PSJ) 

PRMSR-PR ( J  J/PRMS 

WRITE (NW, 193) J, PR ( J),PRMSR 

21  CONTINUE 

DO  22  J ■ 1  *  NX  M 
DO  22  I«1,NTM 

22  P(J, I ) «P  C  J • I )— PA ( J I 
NYM*  8000 

ANN-NYM 

NKl-NXM-1 

DO  24  KK«1,NKL  . 

DO  25  K-1,500 

Kl-K-1 

SUM-0. 

SM-O.SSKKM-O. $SKK0MK-0 .SSKKMK-O. 
SUM1-0. 

DO  40  M«1,NYM 

SUM-SUM  +  P(KK,M)*P (KK»M+K1  ) 
SUM1»SUM1+P( 1  * M) *P  ( KK  +  lf  M+Kl ) 
SKKM«SKKM+P(KK,M)**2 
SKKMK«SKKMK  +  P(KK,M+K1  )**Z 
S  M  «  S  M  ♦  P  { 1  >  M ) 4*2 
SKK0MK-$KK0MK+P(KK+1,M+K 1)442 
40  CONTINUE 

SKKM-SORT  (S*<KM/ANN) 

SKKMK»SORT( SKKMK/ANN) 
SM«SOPT(SM/ ANN) 


SKKQMK-SORT ( SKKOMK/ANN ) 
SUMX(K)»SUM/ANN/SKKM/SKKMK 
SU^Y (K ) «SUM1/ ANN/SKKOMK/SM 
25  CONTINUE 

WRITE  INW,  191 ) 

WRITE(NW>200) (SUMX(K)>K*1>  50  0) 
WRITE(NU,191) 

WRITE ( NW, 200) (SUmY(K)»K*1»500) 

24  CONTINUE 

191  FORM AT ( 1H1 ) 

192  FORM AT (//>5X>" UNCORRECTED  PAVGC 

193  FORMAT  (10X, "CORRECTED  PRMSCSI3, 
200  FORMAT (5F20. 6) 

JXY«5500 
DO  93  J»1>NXM 

WRITE (1) ( P ( J  »  KJ )>KJ*1>16500) 

93  CONTINUE 
1000  STOP 
END 


1 3>  " )  IS'SE20.6) 

)  IS'S  E  20. 6> "  RATIO'S  F10.3) 


ism 
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